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We address stationary patterns in exciton-polariton condensates supported by a narrow external 
pump beam, and we discover that even in the absence of trapping potentials, such condensates may 
support stable localized stationary dissipative solutions (quasi-compactons) , whose field decays faster 
than exponentially or even vanishes everywhere outside the pump spot. More general conditions lead 
to dissipative solitons which may display bistability. The bistability in exciton-polariton condensates, 
which manifests itself in simultaneous existence of two stable and one unstable localized solitons 
with different amplitudes, widths, and exciton-photon fractions under the same physical conditions, 
strongly depends on the width of pump beam and is found to disappear for sufficiently narrow pump 
beams. 

PACS numbers: 67.10.Jn, 03.75.Lm, 71.36.-l-c, 05.45.Yv 



I. INTRODUCTION 

Triggered by a series of remarkable experiments^"— , 
during the last decade one observes a rapidly grow- 
ing interest in study of highly populated (macroscopic) 
exciton-polariton states, which are also identified as 
quasi-particlc Bosc-Einstein condensates (BECs)^'^. Un- 
like the atomic BECs, the condensates emerging from 
the interaction of light and matter are nonconservative 
and exist due to the intense energy exchange with laser 
pump beam, balanced by relatively strong losses of ex- 
citons and photons arising, for example, due to material 
and cavity imperfections. The interplay between gain 
and losses makes such systems an ideal tcstbed for many 
fundamental phenomena in open quantum systems and 
in nonlinear science (in addition to the well appreciated 
practical relevance of physics of microcavities^) . 

The formation of nonlinear dissipative patterns in 
semiconductor microcavities may be accompanied by the 
bistability— that arises due to Kerr-like nonlinearity in- 
duced by the exciton-exciton interactions and that was 
observed for extended pump beams in various setups^"—. 
The theory of this effect was developed within the frame- 
work of different meanfield models based either on equa- 
tions for the polariton occupation number—, where the 
fractions of excitons and cavity photons are connected 
through the Hopfield coefficients^^, or on coupled equa- 
tions for excitons and photonsi^. Bistability of coherent 
spin ensembles described in the two-mode approxima- 
tion by the spinor Gross-Pitaevskii equation with exter- 
nal pump has been reported too^i. Under appropriate 
conditions bright solitons may form in semiconductor mi- 
crocavities as predicted in RefJi and observed in Ref J^. 
All such solitons, reported so far, were detected only on 
a broad constant background due to the use of very ex- 
tended pump beams, whose width is comparable with 
the cavity size — or considerably exceeds a width of lo- 



calized excitation propagating on a broad background — . 
However, in experiment the pump beam may be inho- 
mogcneous and strongly localized in space what leads to 
highly nontrivial collective dynamics of polaritons within 
pump spot — and can be used for observation of rich dy- 
namical phenomena, like direct visualization of quantum 
field patternsi^ or switching in polaritonic transistors^^. 



The possibility of formation of strongly localized non- 
linear excitations in condensates supported by a narrow 
pump beam and the impact of the width of pump beam 
on the localization degree of stationary states have not 
been explored, so far. In particular, the existence of 
highly localized dissipative quasi-compactons, i.e. sta- 
tionary compact waves which decay faster than expo- 
nentially or in some approximation are even are exactly 
equal to zero everywhere except the pump spot, remains 
an open issue. Notice, that the concept of compacton was 
first introduced for continuous systems obeying nonlinear 
dispersion^. Later on compactons have been found in 
discrete — and continuous^ nonlinear conservative lat- 
tices. To the best of our knowledge, dissipative com- 
pactons (representing stable attractors excitable from a 
broad range of initial conditions) were not detected be- 
fore not only in polariton condensates, but in any dissi- 
pative system. 

This work is devoted to the investigation of dissipa- 
tive solitons in microcavities with localized pump. We 
detected the possibility of coexistence of several stable 
solitons with various fractions of excitons and photons in 
the condensate for a certain range of frequencies of pump 
beams, i.e. the existence of bistability, which, however, 
completely vanishes when the width of the pump beam 
becomes sufficiently small. Moreover, for the first time 
to our knowledge, we obtained strongly localized quasi- 
compactons in dissipative exciton-polariton condensates. 
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II. THE MODEL 

We restrict the consideration to a quasi-one dimen- 
sional condensate of strongly coupled excitons and pho- 
tons (designated by the sub-indexes x and c respec- 
tively), which can be described by the Hamiltonian^ii^ 
H = + He + Hni + Hint + Hp where H^^^ are 
the Hamiltonians of the excitons and cavity photons 
given by H^,^c = / dx^^l {fvuj^^c - (tf /2m^^c)d'^ /dx"^) 4-^, 
Hni = (%/2) / dx^'J.^'J.^'a-^'a;, describes the interaction 
between excitons characterized by the strength g > 0, 
Hint = fi^R J dx '^I'^c + h.c. describes the interaction 
between light and matter which is quantified by the Rabi 
frequency 0_r, and Hp — hj dxF{x)e''"''*'^l + h.c. de- 
scribes pump of the photons by an external source, which 
is considered monochromatic and normal to the surface, 
i.e. F{x) is considered real. 

We adopt the meanfield description where the field op- 
erators *x,c are replaced by the respective order parame- 
ters which we represent in the form <j)x^c{x,t)e^'^''* . Then, 
the system of the meanfield equations can be written in 
the form^ii^i2^ 

^lx(Px + ^l(Pc + \(Px\ (Px (la) 



dr 



dr 
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Here we introduced the dimensionless independent vari- 
ables r = \ flc\t and rj = y^md^lcl/hx, where flc = i-^p—^^c 
is the pump frequency detuning from the photon fre- 
quency at fc|| = and we have limited the consider- 
ation to the negative detuning ilc < 0. We also de- 
fined the normalized pump detuning from the exciton fre- 
quency e = {ujp — ujx)/\^c\, dimensionless Rabi frequency 
n = J7/?/|r2c|, and the relation between the masses of 
cavity photons and excitons k = mc/rrix ~ 10~^. In ([T]) 
the normalized cavity losses of photons and the exciton 
decay rates are determined respectively by jc and jx and 
the pump shape is now desribed by /(?;) = F{x)/h\nc\- 
The dynamics of the condensate is determined by sev- 
eral experimentally controllable parameters. Below as 
control parameters we use the amplitude a and width w 
of the Gaussian pump beam /(tj) = aexp (^—ri^/w^), as 
well as the normalized detuning of the pump frequency 
ujp from the exciton resonance, i.e. e (experimentally e 
can be varied by applying mechanical stress to the cavity 
that affects exciton states without changing the energy 
of photons^) . 

Due to smallness of the cavity photon mass, and thus of 
the coefficient k, the kinetic energy of excitons is usually 
neglectedii. While we do not use this approximation in 
numerical simulations, it is useful for qualitative discus- 
sion of the results. Indeed, by neglecting the kinetic en- 
ergy of excitons in (fTa|) one obtains the relations between 
the exciton and photon linear densities rix^c = IV'x.cP and 
the phase difference 9 ^ 9c — Ox'- 
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These relations imply that for e > y/2>^x one can find 
several solutions for exciton densities for a given photon 
density, while for e < \/i^x for each photon density there 
exists a unique density of excitons. Eqs. ^ identify an 
important (real) physical parameters 5 and (p, which we 
define by Se^"^ ~ e ~ ij^- Indeed, far from the pump 
spot (formally in the limit \r]\ — > oo) where n^.c — ^ we 
obtain from ^ that the relation between the densities 
of excitons and photons is determined by S: nx/ric — > ^, 
while the phase difference is fixed by — >■ ip. 

For inhomogeneous polariton condensates supported 
by a localized pump it is convenient to pass from the 
local densities rix^c to the total number of excitons and 
photons Ux,c = J \4>x,c\'^dx. In Fig.[T]we show the numeri- 
cally calculated number of excitons in stationary solution 
as a function of different parameters of the system. While 
for relatively wide Gaussian pump beams in certain pa- 
rameter range wc do observe several soliton solutions for 
the same value of control parameter [Fig. [TJa)], the first 
important result of this work is that only one solution re- 
mains when the width of the pump beam decreases. This 
phenomenon is observed in the dependencies of the num- 
ber of excitons on the pump beam amplitude [Fig. [IJa)] 
and on the detuning e [Fig. [TJb)]. Thus, for wide pump 
beams the dependence Ux{a) is "S" -shaped, while for nar- 
row beams Ux(a) is a monotonically growing function. 
Linear stability analysis of the obtained solutions pre- 
dicts that only upper and lower branches of solitons in 
Ux{a) dependence may be stable, while middle branch 
where dUx/da < is always unstable (stable branches 
are shown black, while unstable branches are shown red 
in Fig. [1]). Therefore, in certain parameter range two 
stable solitons coexist with one unstable soliton. 

While the bistability in the dependence on pump beam 
amplitude a resembles, when it exists, the behavior of a 
homogeneous condcnsat o^^'^^ , in our setting it can also 
appear in the dependence on detuning g^l. Indeed, as one 
can observe in Fig. [Ub), the exciton density is not any- 
more a monotonic function of this control parameter, but 
in the bistable region Ux grows (decays) with e on the up- 
per (lower) branch. Here it is relevant to address the real 
physical scales at which the bistable behavior disappears. 
Assuming the frequency detuning Mlc ^ 1 meV and ob- 
serving that the change in the behavior occurs between 
w ~ 2 and w = 3 wc estimate that bistability disappears 
for pump spot sizes between 3 and 5 nm. The threshold 
spot size can obviously be increased by increase of the 
detuning. 

The phenomenon of bistability is also obvious in the 
dependence Uxi'jx) [Fig- He)] that indicates that upper 
and middle branches merge when losses reach critical 
level and only one stable lower branch remains (the de- 
pendence on 7c is qualitatively similar). 

In Fig. [Ijd) we present the ratio of the number of exci- 
tons and photons in condensate as a function of detuning 
e. Comparing with Fig. [Ub) we observe that the bista- 
bility of localized polaritons is related to different exci- 
ton and photon fractions Ux and Uc. In order to clarify 
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the relevance of this observation we recall the fact that 
the exciton-polariton condensate is a mixed state of exci- 
tons and photons whose relative fractions determine the 
condensate properties. Affecting transmitted or reflected 
light the relative fraction of the quasi-particles is also 
relevant for interpretation of the experimental measure- 
ment^. Since the seminal work of Hopficldi^, practically 
for all estimates involving these fractions the Hopficld co- 
efficients are used. In notations of the present work, for 
the relative fractions of excitons \X\'^ and photons |Cp 
of the homogeneous condensate one has 
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This formula docs not display bistability [for compari- 
son with numerically obtained value Ux/Uc in Fig. [Ijd) 
we show the relative fraction jXp/jCp determined by 
dU]. Since the Hopfield coefficients were used in previous 
studies of bistability in homogeneous exciton-polaritons 
BECs (sec c.giii^), as well as in the justification of the 
meanficld models written for the polariton order param- 
eter (like in the models used e.g. in^^), obviously the 
respective meanfield models cannot capture the effects 
related to variations of the relative fractions of the quasi- 
particles. Remarkably, Fig.[ljd) shows that the branch of 
solutions at negative e values is very well reproduced by 
formula ([3]), while strong qualitative differences appear 
at larger detunings: there the exciton fraction grows sig- 
nificantly stronger than this is predicted by the Hopficld 
coefficients. 

Bistable behavior of exciton and photon fractions be- 
comes particularly evident from typical density distribu- 
tions in stationary solutions shown in Fig. [2] for three 
branches of the bistable curve from Fig. [TJa) (see points 
marked by circles). The figure reveals several interesting 
features. First, the fraction of photons is larger (smaller) 
than the fraction of excitons on the upper (lower) branch 
in Fig.[T^ [notice that the lower (upper) branches in pan- 
els a, b and c of Fig.[l]correspond to upper (lower) branch 
on the panel d]. Despite significant differences in frac- 
tions of photons and excitons, both components occupy 
nearly the same spatial domain and the widths of both 
components exceed the characteristic size of the pump 
beam. 

Moreover, the soliton from lower branch (panel a) is 
substantially narrower than soliton from upper branch 
(panel c) , what is reflected in the respective spatial spec- 
tra sown in Fig. [3] Notice the appearance of oscillations 
on shapes and spectra of solitons from middle, i.e. un- 
stable, branch (panel b). Existence of three spectral 
maxima in Fig. |3] indicate on the possibility of resonant 
mode interactions, which can be triggered by any small 
perturbation, what is likely to be the physical reason for 
the instability of the modes from the middle branch. 

The existence of more than one stable state in a system 
with experimentally adjustable parameters implies the 
possibility of hysteresis. Previous theoretical considera- 
tion^i^i^ dealt with hysteresis of condensates supported 
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FIG. 1: Total number of excitons vs (a) the amplitude of 
pump beam at e = — 1 and 'yx,c = 0.1 (the points correspond 
to solitons shown in Fig. [2] below) , (b) the detuning e at a = 
0.7 and 'yx.,c ~ 0.1 (the curves from top to bottom correspond 
to pump widths w = 5, 3, 2, 1), (c) the exciton losses 7a; at e = 
— 1, w = 5, and 7c = 0.1. (d) The ratio Ux/Uc vs frequency 
detuning at a = 0.7, w — 5, and 7c, a; = 0.1 corresponding 
to soliton solution (curve 1) and the ratio Ux/Uc calculated 
using Hopfield coefficients (curve 2). In all cases we used 
Q, = 2. The red segments correspond to unstable branches of 
solutions, while black segments to stable branches. 



by nearly homogeneous pump (what corresponds to suf- 
ficiently wide pump spots in the experiment). In Fig. |3] 
we show the hysteresis of a localized dissipative soliton in 
the polariton condensate. In particular. Fig. Ufa) shows 
jump from lower to upper soliton branch stimulated by 
small increase of the amplitude of pump beam beyond 
the amplitude at which lower branch ceases to exist (the 
input soliton corresponds to right outermost point in 
Ux{a) dependence). Analogously, Fig. lUJb) shows jump 
from upper to lower branch when pump beam amplitude 
was decreased below the threshold for existence of upper 
branch. 



III. QUASI-COMPACTONS 

Localized soliton solutions presented here are dissipa- 
tive solitons with decaying density distributions, which 
are supported by the localized pump beam. Somewhat 
similar excitations may be obtained in nonlinear op- 
tics^. However, the drastic difference between the case of 
exciton-polariton condensate considered here and optical 
settings is that in the case of condensate the pump is ex- 
ternal, rather than parametric. In this case the decay law 
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FIG. 2; The density distributions at a — 0.7, w = 5, e = — 1 
and 7a;, c = 0.1 for lower (a), middle (b), and upper (c) soliton 
branches corresponding to circles in Fig. [ija). 



of soliton outside the pump domain may be completely 
different from the decay law in the free space what opens 
the possibilities for formation of compactons (or quasi- 
compactons), i.e. solutions which are zero (or displaying 
faster than exponential decay) outside the spot domain. 
In order to explain this phenomenon, we again turn to 
the case k = 0. From Eq. © it follows that at 77 — > oo 
the phase and amplitude of the exciton and photon com- 
ponents arc uniquely related by the parameters S and ip. 
Hence, at ?/ — > 00 there exist only two free parameters 
defining the solution. Indeed, by neglecting the nonlinear 
term (which decays much faster than the linear ones) in 
Eq. ([1} with K = one can obtain the asymptotics of 
the solution at 77 — > 00: 



^^^^^smH{^k-^^)\r,-x\)^^ (4) 
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where 
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and ^ = {'-fc + il'^'jx/S^)/\k\ is the decay exponent. The 
constant C is complex. Its amplitude and phase are two 
constants parametrizing the solution. In order to find the 
particular value of C one has to employ the additional 
constrains (see e.gi^). For example, for even solitons 
4>x,c{'n) = 4'x.c{—Ti) the requirement of zero derivative of 
4'x,c{'ii) at 77 = gives two conditions for definition of 
complex constant C via the system parameters. If /(??) 
is nonzero and decays faster than exponentially (as in 
the present work) two situations are possible: if system 
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FIG. 3: Spatial spectra ^-^^^i^k) = ^'^^^^,^(x)e-^^''dx 
corresponding to the profiles shown in Fig. [2] 




FIG. 4: Hysteresis of the dissipative solitons. (a) Switching 
from the lower to upper branch at a = 0.839 (slightly exceed- 
ing the amplitude at which the lower branch ceases to exist), 
(b) Switching from the upper to lower branch at a = 0.613 
(that is slightly below the pump amplitude at which the up- 
per branch ceases to exist). In both cases w = 5, e = — 1, 
7i,c = 0.1, and Q. = 2. Only the modulus of exciton compo- 
nent is shown. 



parameters are such that C ^ 0, the decay of solution is 
exponential, but if C = the solution is still nonzero, but 
it decays in accordance with the law determined by the 
shape of the pump beam /(ry), i.e. faster than exponen- 
tially. Such solutions are dissipative quasi-compactons 
and their shapes obtained from Eqs. ([1]) are presented in 
Fig. [5] for the case of k = 10~^. Quasi-compactons ap- 
pear only on the upper and middle soliton branches in 
the region of positive detuning, when e > (i.e. when 
cjp > uJx)- Typically, the larger the value of e the faster 
is the decay of tails of solution. From Fig. [Ub) one can 
see that upper soliton branches, hence quasi-compacton 
solutions, can be found in the region with positive de- 
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FIG. 5: The wavefunction modulus for exciton and photon 
components for soliton solution from the upper branch at e = 
(a) and e = 1 (b) for a = 0.9, w = 5, ^x,c = 0.1, and $1 = 2. 

tuning only for sufficiently large pump width w. Quasi- 
compactons can be stable. 

Remarkably, not only quasi-compacton, but also ex- 
act dissipative compacton solutions can be obtained by 
choosing proper pump beam shapes in the limit k = 0. 
To illustrate this, let us withdraw the requirement for the 
pump f{ri) to be real and choose 

fiv) - -foSe-''^ + /o[(3 - ilc)Se~"^ + n^] cos^ 

+1^2/3 cos^ 7] [15 - (19 - ijc) cos^ 7]] , (6) 

for \ri\ < it/2 and / = otherwise (/o is a free param- 
eter, also recall that (5e"^ = e — i'^x)- Then, the exact 
compacton solution of Eq. ([1]) with k = reads 

(t)x -fo^ cos^ 77, (j)c = /o cos^ 1] {f^n'^ cos'^ 1] - Se^'^f) 

(7) 

for 1 77! < 7r/2 and <j>x = <Pc = otherwise. We confirmed 
numerically that this compacton solution is structurally 
(i.e. at K = 10~^) and dynamically stable. We also no- 
tice that the suggested solution is not unique one of such 
kinds. As a matter of fact a large diversity of the dissi- 
pative modes with compact support (i.e. equal exactly 
zero beyond a give region) can be constructed, although 
by using even more sophisticated gain profiles. 

Explicit analytical solution of a dissipative compaction 
allows to address one more important issue, which is the 
dissipative nature of the obtained solutions. To this end 
let us first set 7^ = 7^ = 0. Then ©, ^ still hold 
valid and represent a Hamiltonian system subjected to 
an external force f{r]). Such a system can be viewed as a 
composed quasi-particle with one component (excitonic) 
having an infinite mass (since k = 0) affected by the in- 
homogeneous force applied to the photonic component. 
This force is balanced by the force of interaction between 
the components (determined by fi) giving origin to the 
stationary state described by ([7]) with 7c = 7^ = 0. It 
turns out, however, that stationary localized stationary 
solutions of Eq. ([T]) can be found only when at least one 
of the parameters ^x,c is nonzero. In this case one can ob- 
tain one or three soliton solutions depending on the par- 
ticular value of pump amplitude, as shown in Fig. [ija). 



When both 7^: = 7c = stationary states are unstable 
and one observes considerable radiation from pump re- 
gion. 



IV. CONCLUSION 

Summarizing, we showed that spatially localized pump 
supports strongly localized dissipative solitons in exciton- 
polariton condensates. The width of the pump beam de- 
termines whether the system possesses or is free from 
bistability. We found that in the bistability regimes the 
relative fractions of excitons and photons may signifi- 
cantly differ from the predictions for the homogeneous 
condensate where they are determined by the Hopfield 
coefficients. 

Although bistability is considered to be a common fea- 
ture for externally driven systems (it occurs not only 
in exciton-polariton condensates, but also in lasers, as 
well as in various systems described by driven nonlinear 
Schrdinger equation [30] ) , in real experiments it is usually 
manifested in the abrupt jump of the norm (or power) 
of the output state with increase of the pump amplitude, 
so that observation of localized states from certain power 
range usually cannot be performed due to their instabil- 
ity. In this respect, the obtained here fact that by simple 
change of the size of the pump spot the states from the 
above mentioned power range can be stabilized is rather 
interesting. 

For the first time to our knowledge we presented dissi- 
pative compactons, or more precisely quasi-compactions, 
i. e. solutions which decay faster than exponentially. 
Thus it is of interest further study of quasi-compactons 
from the point of view of driving such modes by laser 
beams, of their interactions, their use for switches and ob- 
servation of the hysteresis phenomenon, as well as many 
other problems typical for the nonlinear physics. 

Finally, it should be stressed that the consideration in 
this work was based on the meanfield model consisting 
of the coupled equations for the excitons and photons, 
thus leaving a series of open questions. Among them 
we mention the extension of the theory for spinor model 
of the condensate; analysis of possible excitations of the 
condensation in excited states and associated energy re- 
laxation^^; extension to two-dimensional geometry; dy- 
namical properties and possibility of managing localized 
excitations, etc. 
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